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Abstract. We study the stability of resonance poles in nN Pn partial wave, particularly the 
Roper resonance, by varying parameters significantly within the EBAC dynamical coupled-channels 
model, keeping a good fit to the empirical amplitude. We find that two Roper poles are stable against 
the variation. However, for higher energies, the number of poles can change depending on how the 
parameters are fitted within error bars. We also developed a model with a bare nucleon which forms 
the physical nucleon by being dressed by the meson-cloud. We still find a good stability of the Roper 
poles. 
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INTRODUCTION 

For extracting A'^* information, first, one needs to construct a reaction model through 
a comprehensive analysis of data. Then, pole positions and vertex form factors are 
extracted from the model with the use of the analytic continuation. Therefore, the A'^* 
information extracted in this manner is inevitably model-dependent. Thus, commonly 
asked questions are how model-dependent are the extracted resonance parameters, and 
how precise do data have to be for stable resonance extraction ? These are the questions 
we address at Excited Baryon Analysis Center (EBAC) at JLab[l], within a dynamical 
coupled-channels model (EBAC-DCC) [2]. We focus on the nN Pu partial wave and the 
stability of its pole positions, particularly those corresponding to the Roper resonance. 
In the region near Roper A''(1440), two poles close to the 7cA threshold were found in 
our recent extraction [3] from the JLMS model [4] (JLMS is one of EBAC-DCC model), 
while only one pole in the similar energy region was reported in some other analyses. We 
examine the stability of this two-pole structure against the following variation, keeping 
a good reproduction of SAID single-energy (SAID-SES) solution [5]. 

• Large variation of the parameters of the meson-baryon and bare A^* parameters of 
the EBAC-DCC model. 

• Inclusion of a bare nucleon state: The analytic structure of this model is different 
from the EBAC-DCC model in the region near the nucleon pole [6], . 



DYNAMICAL COUPLED-CHANNELS MODELS 



Here, we briefly describe the EBAC-DCC model and the bare nucleon model. The 
EBAC-DCC model contains kN, r\N and tztzN channels and the tzkN channel has 
;rA, pN and oN components. These meson-baryon (MB) channels are connected with 
each other by meson-baryon interactions {vmb,m'b')', or excited to bare A^^* states by 
vertex interactions (Ymb-i^n*)- With these interactions. The partial- wave amplitude for 
the M(k) +B{—k) M'(k') +B'{—k') reaction can be conveniently decomposed into 
two parts as TMB,M'B'{k,k' ,E) = tMB,M'B^{k,k',E) + t^^^,g,{k,k',E). The first term is 
obtained by solving the coupled-channels Lippmann-Schwinger equation with vmb,m'b' 
only. The second term is associated with the bare A^* states, and given by 

^MS,M'5'(^,^'>^) = Y.^MB^N*{k,E)[D{E)]iJm^M'B'{k' ,E), (1) 

where the dressed vertex function TN*.^M'B'{k,E) is calculated by convoluting the bare 
vertex Fyv^^M's'l^) with the amplitudes tMB,M'B'{k,k' fE). The inverse of the propagator 
of dressed A'* states in Eq. (1) is 

[D-\E)]ij = iE-m%.)5ij-i:ijiE), (2) 
where m^* is the bare mass of the i-th N* state, and the N* self-energy is defined by 

= q^dqT]^*^MB{q,E)GMB{q,E)TMB^NT{q,E) , (3) 

MB-'^MB 

where Gmb is the meson-baryon propagator, and Cmb is the integration contour in the 
complex-^ plane used for the channel MB. 

To examine further the model dependence of resonance extractions, it is useful to 
also perform analysis using models with a bare nucleon, as developed in Ref. [8]. 
Within the formulation of EBAC-DCC model, such a model can be obtained by adding 
a bare nucleon (A^o) state with mass and A'o MB vertices and removing the direct 
MB N M'B' in the meson-baryon interactions vmb.m'b'- All numerical procedures 
for this model are identical to that used for the EBAC-DCC model, except that the 
resulting amplitude must satisfy the nucleon pole condition: 

E-m'^-L{mM) 

with 

niN = m% + t{mN) and F^mikon) ^ F^^^'{kon) ■ (5) 

Here we have used the on-shell momentum defined by £■ = y^m^f-K^ + aA^I^K^- 

Also, £(mA^) is the self-energy for the nucleon. More details for the calculational proce- 
dure following Afnan and Pearce is found in Refs. [1, 8]. 



RESULTS 



Now we show our numerical results to examine the stability of the Pu poles. We present 
results from various fits by varying the dynamical content of the EBAC-DCC model, 
and by using a model with a bare nucleon. We show in figures the quality of fits of these 
models, and in Table 1 the pole positions from the models as well as per data point 
(Xpd)- fi^*^ poles with the method of analytic continuation discussed in detail in 
Refs. [3, 9]. In Table 1, we also present pole positions from JLMS[4] and SAID-EDS 
(energy-dependent)[5] . 

First we varied both the parameters for the meson-baryon interactions (vmbm'b') 
parameters associated with bare A'^* states (m^*, ri\f*^MB) within EBAC-DCC model. 
The obtained meson-baryon interactions are quite different from those of JLMS. We 
obtained several fits which are different in how the oscillatory behavior of SAID-SES 
amplitude for higher W is fitted. The results from the 2A^*-3p (dotted curves) and 2N*-4p 
(dashed curves) fits are compared with the JLMS fit (solid curves) in Fig. 1. The resulting 
resonance poles are listed in the 3th and 4th rows of Table 1. Here we see the first two 
poles near the nA threshold from both fits agree well with the poles from the JLMS fit. 
This seems to further support the conjecture that these two poles are mainly sensitive 
to the data below W ^ 1.5 GeV where the SAID-SES has rather small errors. However, 
the 2A^*-4p fit has one more pole at Mr = 1630 — i45 MeV. This is perhaps related to 
its oscillating structure near W ~ 1.6 GeV (dashed curves), as shown in the Figs. 1. On 
the other hand, this resonance pole could be fictitious since the fit 2N*-3p (dotted curve) 
with only three poles are equally acceptable within the fluctuating experimental errors. 
Our result suggests that it is important to have more accurate data in the high W region 
for a high precision resonance extraction. 

Next, we show our results obtained with the bare nucleon model, and then address 
the question whether difference in the analytic structure of the TtN amplitude below kN 
threshold strongly affects the resonance extractions. The bare nucleon model is fitted to 
SAID-SES, and at the same time, to the nucleon pole conditions Eq. (5). Meanwhile, the 
original EBAC-DCC model has different singular structure below the ^rA'' threshold. The 
question is whether such differences can lead to very different resonance poles. Our fit 

TABLE 1. The resonance pole positions Mr for Pu [listed as (ReMR, — ImMR) in 
the unit of MeV] extracted from various parameter sets. The location of the pole is 

specified by, e.g., (s7rAf,Sj)A?,S7r;iAf,.S;rAi*p'Vj*o"iv) ~ {upuupp), where p and u denote 
the physical and unphysical sheets for a given reaction channel, respectively, x^d x} 
per data point. 
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(1364, 105) 
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FIGURE 1. The real (left) and imaginary (right) parts of the on-shell Pn amplitudes as a function of 
the 7tN invariant mass W (MeV). A is unitless in the convention of Ref. [5]. 



of the bare nucleon model agree very well with JLMS below W = 1.5 GeV, while their 
differences are significant in the high W region. The corresponding resonance poles are 
given in Table 1 . We also see here that the first two poles near the nA threshold are close 
to those of JLMS. Our results seem to indicate that these two poles are rather insensitive 
to the analytic structure of the amplitude in the region below nN threshold, and are 
mainly determined by the data in the region niN + m^r < W < 1 .6 GeV. 
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